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In this numerical study on two-dimensional Rayleigh-Benard convection we consider 10^ < Ra < 
10^^ in aspect ratio 0.23 < F < 13 samples. We focus on several cases. First we consider small 
aspect ratio cells, where at high Ra number we find a sharp transition from a low Ra number branch 
towards a high Ra number branch, due to changes in the flow structure. Subsequently, we show 
that the influence of the aspect ratio on the heat transport decreases with increasing aspect ratio. 



CN| although even at very large aspect ratio of F 10 variations up to 2.5% in the heat transport as a 

^-H function of F are observed. Finally, we observe long-lived transients up to at least Ra = 10^, as in 

certain aspect ratio cells we observe different flow states that are stable for thousands of turnover 

times. 

^ PACS numbers: 47.27.-i, 47.27.te 
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FIG. 1. a) Nu versus Ra for Pr = 1 in aspect ratio F = 0.23 (orange circles), F = 1/3 (lime diamonds), F = 1/2 (green 
squares), F = 2/3 (blue upward pointing triangles), F = 1 (purple pentagrams) and F = 2 (red downward pointing triangles) 
samples, b) Re2 versus Ra for Pr = 1 and in F < 1 samples. The averaging time ranges from 4000 t^; to 200 ts- 



I. INTRODUCTION 

In Rayleigh-Benard (RB) convection [1, 2 a fluid in a box is heated from below and cooled from above. This 
system is paradigmatic for turbulent heat transfer, with many applications in atmospheric and environmental physics, 
astrophysics, and process technology. Its dynamics are characterized by the Rayleigh number Ra = /3gAL^ / (kv), the 
Prandtl number Pr = u/hi and the aspect-ratio T = d/L. Here, L is the height of the sample and d its width, /3 is 
the thermal expansion coefficient, g the gravitational acceleration, A the temperature difference between the bottom 
and the top of the sample, and v and k the kinematic viscosity and the thermal diffusivity, respectively. 

To better understand geo- and astrophysical phenomena it is important to know the heat transfer in the high Ra 
number regime. Nowadays, most experimental and numerical measurements of the heat transfer, indicated by the 
Nusselt number Nu^ in 3D agree up to Ra ^ 2 x 10^^ (see Ref. [1 for detailed references) and these data are well 
described by the Grossmann-Lohse (GL) theory [3]-[6]. Unfortunately, there are considerable unexplained differences 
in the heat transfer measurements for higher Ra [T-Tl . 

Explaining the differences between the high Ra number measurements is subject of an ongoing discussion [TJ[TTl[T4l 
[201 l2T] . The difficulty is that experiments in the high Ra number regime are very challenging as it requires very large 
setups that operate under high pressures [iS^-TlQ^, or use liquid helium near its critical point [T Him [T2UT4] . Therefore 
only a few setups in the world can perform high Ra number experiments and testing the influence of some effects, e.g. 
the influence of the aspect ratio, is very time consuming [M]. In 3D simulations it is difficult to study the problem, 
as these simulations are very CPU-time intensive [2T] . 

One possibility to limit the required computational time is to study the influence of several effects in 2D RB 
convection in order to select the most interesting cases that need to be studied in 3D. One may wonder whether 
the dynamics of 3D RB convection are sufficiently captured in 2D simulations. This question has been addressed by 
Schmalzl et al. [22] and they showed that in relatively large aspect ratio cells, where several convection rolls can form 
next to each other, and for Pr > 1 2D RB convection global quantities such as the Nu number and the Re number 
show similar behavior in 3D. Therefore 2D simulations have been used to investigate a number of aspects observed in 
3D convection [HI ESHSl |30l |3T] . 

It has been argued that the differences between high Ra number experiments are due to the formation of different 
turbulent states. Experimental indications for coexistence of different turbulent states in RB convection are found 
by several authors [TTJ [131 [13 [HI I32H341 [43] and are discussed theoretically by Grossmann and Lohse [36] . It has 
been shown by Xi & Xia [34] and Weiss & Ahlers [43] that in a F = 0.5 sample the flow can be either in a single 
roll state (SRS) or in a double roll state (DRS), each with a specific heat transport. Later van der Poel et al. [37] 
confirmed this effect in 2D simulations and showed that the flow state, and thus Nu^ can strongly depend on the 
aspect ratio. Experiments by Roche et al. [M^ show that the transition towards a regime in which a strong heat 
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FIG. 2. Rttc versus F, both on logarithmic scale. A least-squares fit, which is indicated by the black solid line, gives Rac ^ F 
The blue and red dashed lines give Rac ^ F~^ and Rac ^ F~^, respectively. 



transfer enhancement as function of Ra is found strongly depends on the aspect ratio and the experiments by Niemela 
& Sreenivasan [11 reveal the presence of a low and a high Ra number branch. 

Using a direct numerical simulation scheme for 10^ < Ra < 10^^ we show that also in small aspect ratio 2D RB 
convection there is a low and a high Ra number branch for the heat transport and that the transition at the critical 
Rayleigh number Rac between the branches increases with decreasing aspect-ratio. Subsequently, we study the heat 
transport as function of the aspect ratio at fixed Ra. Bailon-Cuba et al [38] have studied the effect on Nu of 
increasing the aspect-ratio in 3D numerics and found that Nu depends on the number of horizontally stacked rolls. 
Now, we try to estimate how large the aspect ratio should be before finite size effects become unimportant. 

The objective of this paper is to study the phenomena associated with varying Ra and F to a large extent. The 
corresponding observations of the double branching in Nu{Ra) and vanishing variations in Nu{T) at large Ra and F 
can be found in chapters III and IV, respectively. Furthermore, in chapter V we will discuss some long-lived transient 
effects we observe at Ra = 10^ and Pr = 0.7. These three cases have overlapping origin and secondary effects, such 
as multi-stability and roll break-up. 

In the case of long-lived transients, the exhibited multiple roll states with different corresponding Nu and Re, 
survive for such long times, that the roll state can either be interpreted as a transient or as stable. This classification 
is very difficult, as the insuperable finite run time of simulations makes it impossible to ascertain that the observed 
state is not a transient. To overcome the ambiguity between transient and multi- stability, the cases for which this 
classification is unclear will be termed as transients and these will be annotated with the simulation run time. This 
time is an indication of the stability of the observed states. We conclude the paper with a discussion of the presented 
results. 



The code on which the results in this paper are based is a fourth order finite difference discretization of the 
incompressible Oberbeck-Boussinesq equations. The velocity boundary conditions on the walls are no-slip. The 
(relative) temperature is fixed at A/2 at the bottom plate and —A/2 at the top plate, with adiabatic sidewall 
conditions. The code has been described and tested in detail by Sugiyama et al. [23-25 and this code was used to 
explain the non-Oberbeck Boussinesq effects observed in 3D experiments of Ahlers et al. [24 , and Zhou et al. j27l[28] 
and Sugiyama et al. [39] found a good agreement between results obtained using this 2D code and experiments 
performed in quasi-2D RB samples. The grid resolution we used obeys the strict criteria formulated in refs. [20) |4Q]. 
The Nusselt number Nu is calculated with four different procedures. The Nu used in this study is the average of 
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FIG. 3. Flow visualisation for the flow in a F = 0.23 sample with Pr = 1 at (a) Ra = 10^, (b) Ra = 10^, (c) Ra = 10^°, 
(d) Ra = 10^^, and (e) Ra = 10^^. The velocity field is represented by the vectors and, red and blue indicate high and low 
temperature, respectively. Note that the colormap is not constant. For Ra =10^^ the activity of the rolls is increased and at 
Ra — 10^^ the roll structure appears to be completely broken down. Flow field movies corresponding to the sharp transition 
can be found in the supplemental material. 



evaluated at the top and bottom plate. Here {•)A,t denotes the average over a horizontal line and time. The resulting 
Nu is checked using the relation between the volume- averaged dissipation rates and Nu [29 , given by 

Nu = ^Ra-'Pr^{idiUif)v,t + 1, (2) 



Nu=^{{dmv^„ (3) 

where u is the velocity normalized with the free fall velocity U = ^/Pg/S.L^ is the dimensionless temperature < ^ < 1 
and {•)v,t denotes the average over the complete domain and time. As the Nusselt number should not be dependent 
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on the procedure by which it is calculated, comparing the results of the different procedures is an indication of the 
numerical error. For all data points in this study the error calculated by comparing the results is smaller than dot 
size. The time convergence of Nu is checked by comparing the time- averages over half the duration of the simulation 
and the full duration. Again, the resulting error is for most cases smaller than dot size. However, this error is not 
necessarily monotonically decreasing with increasing simulation time due to slow trends and rare events. For an 
indication of the uncertainty associated with the calculation of the numerical quantities, the typical averaging time 
will be annotated with the data. The used unit of time is the large eddy turnover time, given by Ie = 47r/(|cJc|)t, 
where Uc is the vorticity at the center of a roll. 



III. HIGH RAYLEIGH NUMBER CONVECTION 



Figure [T^ shows the heat transport as function of Ra for different aspect ratio cells. Just as in the experiments of 
Niemela & Sreenivasan [llj a low and a high Ra number branch is observed and the Nu difference is larger for the 
low aspect ratio cells. One can see that the Ra number at which the transition from the low to the high branch sets 
in increases with decreasing aspect ratio. To find the transitional Ra number Rac we extrapolated the low branch, 
which scales approximately as Nu ~ i?a^'^^, and fitted the points found in the transition region with a power law. 
The intersection of the resulting lines gives Rac^ which is calculated for 0.23 < F < 2/3. Figure |2] shows that Rac 
scales approximately as Rac ~ ^-2-44^ although the F-range is of course very limited. This is very similar to the 
result obtained by Roche et al. p!4l, who performed high Ra number experiments in different aspect ratio samples 
and observed a sharp increase in the heat transfer, which sets in beyond a critical Ra that scales as Rac ^ r~^-^^. 
However, we note that both data sets would also be in reasonable agreement with a F~^ scaling, which would be 
equivalent to saying that a Ra number based on the diameter d and not the height L of the cell is constant p!4] . 

The transitional Ra numbers we find are well below the value at which the transition towards the ultimate state 
is expected [l] and therefore it is important to know what effect is responsible for the sharp increase in the heat 
transport. To answer this question we visualized the flow field for the simulations in the aspect ratio F = 0.23 
sample, where the transition is most pronounced. Figure |3] reveals that for low Ra there are vertically stacked stable 
rolls. The flow structure is substantially different in the higher Ra number cases. For Ra = 1 x 10^^ the increased 
temperature gradients in the corner rolls indicate that they are more active than at lower Ra^ i.e. the corner rolls 
exchange more thermal energy with the large scale flow. As this increases the effective plate area that contributes to 
the heat transport, this leads to heat transport enhancement. For even higher Ra the rolls (partially) break up due 
to increased buoyancy forces. The plumes break up the rolls and at the same time have a decreased travel distance 
towards the opposite plate. This leads to a higher heat transport. In order to verify this picture we calculated Re2^ 
where Rci = uf-^^L/u is the Reynolds number based on the root mean square velocity with the index i = 1,2, a 
indicating horizontal, vertical and absolute velocity, respectively. Figure [l] shows that the transition in the Nu number 
coincides with a strong increase in the vertical velocities. Therefore we can conclude that the transition we see here 
is due to a transition between different turbulent states - the breakdown of the corner-roll - and is not a transition 
towards the ultimate regime. 

When we compare the 2D results with the measurements of Roche et al. [14] we notice that the Ra number at 
which the sharp increase in the heat transfer is observed happens at a lower Ra in 2D than in 3D. Another difference 
is that we observe a low and a high Ra number branch for the heat transfer in 2D, whereas no high Ra number branch 
is observed in the experiments of Roche et al. [14]. However, we note that Niemela & Sreenivasan [11] did observe 
two branches and it could well be that this second branch could also be found in the Roche et al. [14] experiments 
at higher Ra. Although one should be careful with comparing the results from 2D and 3D RB convection in small 
aspect ratio samples, as we will discuss in detail below, we emphasize that it is important to keep in mind is that not 
every increase in the scaling of the heat transport as function of Ra that is observed indicates a transition towards 
the ultimate regime, as we showed here for the 2D case. 

Again we stress that one should be careful with a direct comparison between 2D and 3D results in small aspect 
ratio cells. For small aspect ratio cells the dynamics in 2D and 3D RB convection can be different as vertically stacked 
rolls are formed. In 2D vertically stacked rolls can lead to a larger reduction of the heat transport than in 3D RB 
convection as warm (cold) fluid in the bottom (top) part of the cell is more likely to be trapped than in 3D. In fact 
experiments and simulations of 3D RB convection, e.g. Bailon-Cuba et al. [S^, show a weaker aspect ratio dependence 
than observed in the 2D simulations (section [Tv| ) . 
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FIG. 4. a) Nu as a function of T for Pr = 4.3 and Ra = 10^. The depicted 1 < T < 12 range captures the SRS up to the 
horizontally stacked lORS (n = 10), each represented by a cluster of points from low to high F. b) n/F approaches a constant 
value of approximately 0.77. c) The function ^(F) = Re2/ (Rein/F). The error on the time-convergence is smaller than the dot 
size as the averaging time is more than 1000 Ie- 



IV. LARGE ASPECT RATIO 



To find the critical horizontal length of the system at which the heat transport becomes nearly independent of 
the aspect ratio we extend the dataset for Ra = 10^ and Pr = 4.3 of van der Poel et al ^ towards large F to 
include horizontally stacked rolls and the corresponding transitions between them. The most remarkable observation 
is that the transitions between the horizontal roll states have a different shape than the transitions between vertical 
roll states, see figure [4^ and Ref. [3T. The amplitudes of the Nu variation over the transitions are surprisingly large, 
especially the transition between the SRS and the horizontal DRS. What cannot be extracted from this figure is that 
the lifetime of the transient states increases close to a jump. 

Figure [4]3 reveals that the ratio /(F) = n/F approaches a constant value of approximately 0.8 at large F, where n is 
the number of large scale rolls in the system. Here we define that for vertically stacked rolls n < 1 and for horizontally 
stacked rolls n > 1. Comparing with Nu{T) in figure [4^, it can be seen that there is a strong correlation between Nu 
and n/F. This indicates that the asymmetric deformation of the rolls is connected to Nu, which can be demonstrated 
using the concept of plumes driving the bulk flow and transporting the heat. When a roll is horizontally elongated, 
the ratio between the number of rolls n and F decreases, this results in a lower Nu as the number of thermal hotspots 
per horizontal plate length decreases as 1/F for fixed n. This explains the similarity of the Nu{T) profiles for fixed 
n between the Nu{T) plot in figure [4^ and the /(F) plot in[4]3. In addition, it appears that the amplitudes of the 
Nu jumps between these plots are similar. The amplitudes decreases for each consecutive transition from n to n + 1. 
This can be explained by using identical reasoning for the profile for fixed n, applying to a transition. We then get 
that the jump amplitude ANu^'^'^ between a n and n + 1 roll state scales as 
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FIG. 5. a) The amplitude of the jumps between different in roll states in Nu{r), ANuJ^^^ and the C/(r + D) scaling as a 
function of n. Here the scaling constants C = 6.9 and D = 1. b) The Reynolds number Rea as a function of F appears to 
asymptotically approach a constant Rea, oo value at high F. The fit with equation ([5| is plotted as the dashed line. 
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FIG. 6. Nu versus Rei (a) and Re2 (b) for Ra = 10^ and Pr = 4.3. Different n are indicated by the colored symbols, see 
legend figure [5b. A negative correlation between Nu and Rei and a positive correlation between Nu and Re2 can be seen. 



AiV<- ~ ^ - f = f • (4) 

In figure ANu'^'^'^^ extracted from Nu{T) in figure [4^, is plotted together with a scaled function C/{T-\-D). The 
offset in the denominator stems from the fact that the reasoning behind the derivation of the scaling is only valid for 
n > 1. The agreement between these plots is very good and therefore the fit can be used to extrapolate ANu]^'^^ to 
higher F. Using the asymptotic values for Nu and n/F, which are 24.92 ± 0.06 and 0.77 ± 0.01 respectively, the Fc 
where the jump magnitude is less than 1% of Nu is found to be Fc = 26 ± 1, which corresponds to n ~ 20. 

Not only Nu{T)^ but also Rea{T) can be used for extrapolating. Figure [sJd shows that Rca increases with F and 
approaches an asymptotic Rca^oo- Similar to NuiV) and for n > 1, the trend is well approximated by 



i?e{^^(F) = i?e,,oo(l-0.01^), 



(5) 
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where Tc is the aspect-ratio where the function Re{l (F) is 99% of Rea^oo- Using the F's just after a jump results 
in Fc = 22 ± 1, which is close to the prediction based on Nu(r), considering the quantity used for extrapolation is 
very different. 

The loss of F dependence is accompanied by an increase in multi-stability of the roll states. For large F the transients 
for F's near a jump in n become very long-lived. As an example, for F = 13.6 both the n = 10 and n = 11 roll states 
survive for longer than 3000 tE- For comparison, a typical transient state for this Ra = 10^ and Pr = 4.3 is less 
than 100 Ie- The corresponding Nu differs substantially between these roll states at Nu = 24.68 and Nu = 25.11 for 
n = 10 and n = 11, respectively. Thus, even though the Nu dependence on F decreases monotonically for larger F 
per roll state, which can be seen in the Nu{T) plot shown in figure [4j at large F different roll states have overlapping 
stable regions in F, meaning that hysteresis effects become important. Thus, in spite of the proper extrapolations, 
multi- stability might prevent a solid indication for a Fc, at which the F dependence of global properties is lost. 

The effect of horizontal roll asymmetry is further highlighted by the plots in figure |6] The deformation indirectly 
results in decreasing Nu by decreasing Re2. Figure (6)3 shows that Nu and Re2 are positively correlated. For F > 1.1, 
Re2 decreases monotonically with increasing F due to roll asymmetry up to the F where n jumps to n + 1. In contrast, 
Nu is negatively correlated with i?ei, see figure |6^, which can be explained by a flux conservation argument similar to 
the one presented by Grossmann and Lohse [41^. Considering the LSC to be of simple elliptical shape, the horizontal 
flux has to be equal to the vertical flux. Normalizing with the height of the system, the flux balance is 

Indeed, figure [4]3 reveals that this function is in the order of 1 over the studied F range, which explains the negative 
correlation between Nu and Rei by relating Rei to Re2. 



V. LONG-LIVED TRANSIENTS 

This chapter differs from the previous two in the sense that it does not vary a single control parameter to an large 
extent. Instead, a phenomenon is presented in the form of very long-lived transients. Van der Poel et al. [37 showed 
that for Ra = 10^ and Pr = 0.7 the final state and thus Nu can depend on the initial conditions of the flow field. 
Whether this feature of moderate Ra number RB convection can be found for high Ra number flows is an important 
and open question. Multiple states have been reported for high Ra 3D experiments by Roche et al. [l^ and Weiss 
& Ahlers [43], where the system dynamically switches between two meta-stable states, each with a different Nu. In 
addition, Ahlers et al. [44 showed multiple Nu{Ra) scalings, dependent on minor external temperature variations. 
However, many details of multiple states in RB flow and their effect on Nu are still unknown. 

We now analyze multi-stable states at Ra = 10^ and Pr = 0.7 in 2D, where the largest roll in the flow is in general 
significantly smaller than the limiting dimension of the enclosing boundaries. This qualification implies that simple 
roll states are broken down. These turbulent states with broken-down rolls are named uncondensed, where the notion 
of condensation is followed from [45 . Here, a state is named condensed when the largest roll in the system is of system- 
size, i.e. is spectrally condensed with the energy piled-up in a scale close to system-size. In RB flow the dynamics of 
the coherent structures differ substantially between condensed and uncondensed states. In the uncondensed case the 
largest rolls are no longer stationary on a time scale shorter than a turnover time. Instead, they are mobile and cause 
thermal plumes to move through high-shear central regions, which is very different from the condensed states where 
plumes crawl along the sidewall. Whether the uncondensation is an effect of the inverse energy cascade or a result 
of a growing buoyancy forces compared to inertial forces or both, is not known and beyond the scope of this paper. 
Nevertheless, the condensation terminology is used in this chapter to highlight the similarity between condensation 
in two-dimensional RB flow and classical 2D [46] and thin layer [45j turbulence. 

The uncondensed roll states for high Ra number flows are composed of small and mobile rolls. Therefore an 
advanced identification algorithm is required to obtain statistics on the roll state for these cases. In the remainder of 
this section we first discuss the vortex detection algorithm of van der Poel et al. [37 in detail before we apply it to a 
set of simulations at Ra = 10^ and Pr = 0.7, where very long-lived transients are found. 



A. The vortex detection algorithm 

In order to analyze the large scale flow structure we look at the second invariant of the velocity gradient tensor 
Aij = djUi. The criterion of a positive second invariant [47l[48], formally known as the Q-criterion, is used to identify 
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(a) (b) (c) 



FIG. 7. Example vortex detection of a flow field snapshot for Ra = 10^, Pr = 4.3 and F = 0.85. (a) Flow visualization with 
the temperature and velocity indicated by color and vectors respectively, (b) Mapped by vortex criterium and sign of vorticity. 
Blue and red indicate clockwise and counter clockwise rotation, respectively. Green indicates hyperbolic field, (c) End results of 
the algorithm where the ellipses are colored by direction of rotation. The DRS exhibited by the flow field is captured accurately 
by the vortex detection algorithm. 

vortex regions. For incompressible two-dimensional flow the criterion results in 

A?i + A12A21 < 0. (7) 

A segmentation algorithm is employed to connect the grid points encompassed by the vortex regions after the grid 
resolution is reduced. The purpose of the grid reduction is to decrease the computational cost of the segmentation. 
Although some information is lost in the process, it is only detrimental for small structures, has little effect on the 
LSC structure and the reduction value can be adjusted. The edge points are used in a non- iterative ellipse-fit to model 
the vortices in the flow as ellipses. In general, this shape is more fitting for larger structures and as detecting these is 
the main purpose of this algorithm, it seems justified to use ellipses. The direct output consists of the six coefficients 
of the quadratic conic section. These coefficients are used to calculate the center coordinates, the tilt angle and the 
radii of the ellipse. It requires five floating points and an additional boolean for the direction of rotation, to store the 
information of a single vortex. The operation of the algorithm is visualized in figure [7[ Using this method we can 
now store all information about the main flow structures for long simulations at Ra = 10^ in only 0.15 GB, whereas 
storing the complete flow field over the same time-interval would require 1000 GB. This opens the possibility to study 
rare events of which the information is usually lost. 



B. Different states in the turbulent regime 

Figure [8] shows Nu as a function of F. In this figure the colored symbols indicate the observed roll state using the 
vortex detection algorithm described above. The uncondensed SRS and uncondensed DRS states are distinct from 
each other by their Nu and the number of equal sized largest vortices. A correlation between the roll state and Nu 
can be distinguished as the points are clustered by Nu and the roll state. The simulation we ran for F = 0.950 snapped 
into a SRS state after 750 t^;, unlike neighboring points, which remained in the uncondensed SRS state for longer 
than 2200 tE- Forcing the system in a condensed SRS by using the steady roll state solution of F = 0.950 reveals 
that there is a possibility of multi-stability around this point as this state remains stable for longer than 1800 t^;, 
see figure [8]3. However, for F's further away from F = 0.950, the SRS does not survive and the system immediately 
becomes uncondensed. This happens for e.g. F = 1.100, where this break down occurs within one turnover time tE- 
Thus there appears to be hysteretic behavior not only over the jumps between the different roll states [37], but also 
over the transition to condensed states. 

The SRS has a higher Nu than the uncondensed SRS. This is unexpected, as flow field movies seem to reveal 
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FIG. 8. (a) Nu(r) for different states for Pr = 0.7 and Ra = 10^. The error on the time- convergence is smaller than the dot 
size as the averaging time is more than 1000 ts- (b) Zoom at multi-stability region 0.90 < F < 1.00. 2D PDF plots of the 
X and z coordinates of the vortex centroids for (c) F = 0.73 (uncondensed DRS), (d) F = 0.945 (uncondensed SRS) and (e) 
F = 0.95 (SRS). In (e) the original algorithm output is convoluted with a gaussian for illustrative purpose. This results in the 
peaks being displayed wider than they are. 



enhanced mixing and heat transport due to the presence of the secondary roll in the uncondensed SRS. However, a 
closer investigation suggests that the thermal plumes are not directly transported to the opposite boundary layer in 
the uncondensed SRS, which is the case for the SRS. For the TRS (triple roll state) and DRS, this effect is opposite. 
Here, the break up of the simple roll states increases Nu^ similar to the effect seen in chapter HL 

In contrast, the Nu{r) dependence within the uncondensed SRS, uncondensed DRS and DRS regions appears to 
be absent. However, large Nu variations can be seen at the TRS region, where due to an increased restriction of 
the flow, the largest scale of the flow approaches the length of the limiting dimension of the system and thereby a 
return of strong aspect-ratio dependence. A collection of PDFs of the x and z coordinates of the vortex centroids for 
specific and relevant roll states can be found in figure [S] These plots are based on the output of the vortex detection 
algorithm. In figure [8]i a clear fingerprint of the 'orbiting rolls' uncondensed SRS can be observed as a large central 
peak with a concentric ring. The latter represents the smaller roll(s) rotating around the large central roll. Distinct 
from this figure is the uncondensed DRS PDF in figure [Sj), which reveals that even different uncondensed states can 
be distinguished. In contrast, the SRS PDF has three sharp peaks, which are shown in figure |8^. The PDFs reveal 
that the mobility of the rolls increases for high Ra uncondensed roll states, as the peaks are more smeared out for 
these cases. 
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VI. DISCUSSION AND CONCLUSION 

To summarize, here we considered three cases of two-dimensional Rayleigh-Benard convection in direct numerical 
simulations. First, we showed that in small aspect ratio 2D RB convection there are two branches for the heat 
transport as function of Ra. A lower branch for relatively low Ra and a higher branch for high Ra. By necessity, 
there is a strong increase of the heat transport in the transition region. This transition takes place for higher Ra and is 
stronger when the aspect ratio is smaller. We show that the heat transport enhancement coincides with the transition 
from a flow state with vertically stacked stable rolls to a regime in which the thermal driving becomes strong enough 
to (partially) destroy the rolls. We discuss the similarities and difference with respect to the experiments of Roche 
et at. [14 and Niemela & Sreenivasan [11 and although one cannot directly compare the 2D and 3D results in this 
case it is clear from this example that one should be careful in identifying an increase of the heat transport scaling as 
function of Ra as the transition towards the ultimate regime. 

Second, we showed that the influence of the aspect ratio on the heat transport decreases with increasing aspect 
ratio, although even at a very large F = 10 variations up to 2.5% are observed due to changes in the flow structure. 
Extrapolating RcaiV) towards high F gives Fc = 22 ± 1, where aspect-ratio dependence of global properties is less 
than 1% of their asymptotic values. Using the jumps in NuiV)^ a value of Fc = 26 ± 1 is found. 

Finally, we observe very long-lived transients / different turbulent states at Ra = 10^ and Pr = 4.3 for certain 
very large aspect ratios and at Ra = 10^ and Pr = 0.7. In the latter regime rare random events cause the system 
to switch between states. The probability of these events occurring, the stability of the resulting state and the effect 
on Nu appear to be very sensitive to the control parameters. In particular, the aspect-ratio F has a large influence 
on these phenomena; varying F reveals a rich phase space for a given Ra and Pr. For these parameters the flow is 
so turbulent that the different states with their rare events can only be identified with the use of an advanced vortex 
detection algorithm, which has revealed the possibility of hysteretic behavior over the transition to condensed states. 
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